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Abstract. Stellar limb-darkening is an important tool for constraining the properties 
of a stellar atmosphere. We present a novel method for relating the fundamental stel- 
lar parameters mass and radius to limb-darkening laws using grids of spherical model 
stellar atmospheres. This method is applied to interferometric observations of the red 
supergiant Betelgeuse, where an unique measure of the stellar mass is determined. 



1. Introduction 

The M2Iab star Betelgeuse is one of the brightest stars in the night, and one of the most 
studied. However, for all that is known about the star, there are still a number of myster- 
ies. One such mystery is the mass of Betelgeuse. As a star without a binary companion, 
there exists no direct metho d to measure i ts mass, a mass estimate is possible only using 
theoretical modeling. Dola n et al. fit stellar evolution models to measured val- 



ues of the stellar radius, effec tive temp erature and mass-loss rate to find a stellar mass 



M = 21 +2 Mq. Furthermore. lLobeT& Dupreel fcOOd) used detailed sprectral synthesis 



to derive the surface gravity of log g = -0.5 for Betelgeuse. Coupled with the radius 
assumed by Dplan et al. (2008), that gravity implies a mass M w 9.5 Mq. Clearly, these 
two estimates do not agree. 

The purpose of this work is to present a new measurement of the mass of Betel- 
geuse, using a novel method b ased on observations of the star's intensity profile from 
narrow H-band interferometry dHaubois et al.ll2009l) . the properties of flux-conserving 



limb-darkening laws, and spherically-symmetric model stellar atmospheres. 



2. Model Stellar Atmospheres 



We compute a gr i d of m odel stellar atmospheres using the S Atlas code developed by 



Lester & NeilsonI (I2008r) . which is a spherical version of the Atlas code developed by 
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Kurucz The code assumes local thermodynamic equilibrium and spherically- 



symmetric hydrostatic equilibrium and radiative transfer. The cod e has been tested 
by comparison to spher ically-symmetric Phoenix ( Hauschildt et ali .1999) and Marcs 
( Gustafsson et alj|2008l ) models. Computed intensity profi les have been used to fit K- 
band interferometric observations ( Neilson & Lester I2OO8I) . agreeing with previous re- 
sults (Wittkowski et al. 2004 ). 

The grid is computed in three dimensions, varying effective temperature, gravity 
and stellar mass. The effective temperature range is 3000 - 8000 K in steps of 200 K, the 
gravity range is from log^ = -1 - -i-3 in cgs units in steps of 0.25, and the mass range 
is M = 2.5 - 20 Mq in steps of 2.5 M©. There are approximately 2000 models in this 
grid, although not all combinations of the fundamental parameters produced converged 
models. Model intensity profiles are computed as a function of wavelength and ju, 
where /u is the cosine of the angle formed by the line-of-sight point on a stellar disk and 
the center of the disk. Each profile is computed for 1000 equally spaced ju-points. 



3. Limb-darkening law 

In this work we consider the limb-darkening law 

/ 3 5^ 

— = \-A-B + -Afi + -B^^, (1) 
2'K 2^4^^ ^ ^ 

where 'H is the Eddington flux. This law has two important properties that differ from 
other types of laws. The first property is that the la w conserves stellar flux by definition. 
The second property is of interest for this work. [Fields et aD (l2003l) found that limb- 



darkening fits to intensity profiles from plane-parallel model atmospheres predicted 
relations that all intersected at a fixed point, /ii, with the same normalized intensity. 
The authors compared these profiles to lim b-darkening relations from microlensing ob- 



servations, the observations did not agree. iNeilson & Lesteti (12011a ) showed that this 



fixed point occurs because the limb-darkening coefficients A, and B, are coiTclated such 
that A oc aB. The slope of the correlation is a function of the ratio 77 = j lii-i) ^fild^| J, 
where J is the mean intensity. This ratio is comparable to the Eddington factor used in 
modelling stell ar atmospheres f = K/J, where K is the second moment of the intensity 
(lMihalaslll978h . In plane-parallel model stellar atmospheres, the Eddington factor is 
approximately constant, and similarly ij is also constant. 

In spherically symmetric model stellar atmospheres, the Eddington factor an d 



T] are not constant. This is discussed in greater detail in INeilson & Lesteri (I2011ah . 



Neilson & Lester|(|2011b') discovered that the variation of rj, hence a, and hence the 



fixed point and intensity at the fixed point is correlated to the amount of extension of 
the stellar atmosphere. Furthermore, atmospheric extension can be shown to be a func- 
tion of the ratio of the stellar radius and mass, R^jM^. Therefore, if one can measure the 
limb-darkening coefficients for this law then one has a measure of the atmospheric ex- 
tension and the ratio of the radius and mass. The relation between the values of /?« IM^ 
and the fixed point jiy and the intensity at the fixed point I{fii)/2'}-( is shown in Fig. [T] 
for our models with effective temperature T^s ^ 3000, 4000, 5000, 6000, and 7000 K. 
The fixed point and intensity are functions of R^jM^, such that 

^il=C^{\ogRJM,f + D^, (2) 
l^ = Ci{logRJM,f+D,, (3) 
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Figure 1. The fixed point jji and normalized narrow H-band fit intensity at the 
fixed point /h(//i)/2'K as a function of atmospheric extension parameterized as 
Rt/M,. Red filled circles correspond to Teff - 3000 K models, green open squares 
4000 K, blue open circles 5000 K, magenta downward pointing triangles 6000 K and 
pale blue upward pointing triangles are 7000 K. 



for some given effective temperatures. 



4. Method 

The apparent existence of the fixed p oint in the chosen Um b-darkening law combined 
with interferometric observations from Haubois et al. I (l2009h suggests a method to mea- 



sure the mass of Betelgeuse. The first step is to determine the stellar radius for Betel- 
geuse. We fit the interferometric observations, shown in Fig.[2]along with the predicted 
visibilities from a stellar atmosphere model, to using narrow H-band intensity profiles 
from the model atmospheres to get a best-fit angular diameter. The angular diameter 
can be combined with the measured distance, d = 197 ± 45 pc (Harper et al. 2008), 
and the bolometric flux Fboi = (111.67 ± 6.49) x 10"^^ W cm"^ dPerrin et al.ll2004l') of 
Betelgeuse to predict the stellar radius, and effective temperature. 
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Figure 2. The observed squared visibilities from lHaubois et af] (l2009h are shown 
as crosses as a function of baseline divided by effective wavelength - 1 .64 yum in 
cycles per arcsecond, along with the predicted squared visibility, shown as a dashed 
line, for a T^ff - 3600 K, logg - 0, and M - 12.5 Mq intensity profile. 



The second step is to fit interferometric observations again, but with the limb- 
darkening law to measure the coefRcients, A, and B. Using the correlations from Eq. 2 
and 3, with the limb-darkening law Eq. 1, we have three equations and three unknowns 
that can be solved to determine Rt/M^. Combining this result with the predicted radius, 
then we can measure the mass of Betelgeuse. It should be noted that the coefRcients, 
A, B are degenerate functions of spherical extension because they are both functions of 
J I^JUdpLKin-l) and 7/(2'K). Therefore, one should not directly compare the values of 
the coefficients to the values of /?*/M*. 



5. Results 



We fit the interferometric observations in the same manner as Haubois et al. 1 (l2009h . 



who fit intensity profiles to the observations in the first two lobes of the visibility curve, 
even though the visibility curve is observed to fourth lobe. The authors noted that the 
visibility curve is affected by asymmetries in the photosphere, so it should be ignored. 
We directly fit intensity profiles from the model atmosphere grid to the first two lobes 
of the interferometric observations, and m easure the angu lar diameter of Betelgeuse to 
he 6 = 44.93 ± 0.15 mas, agreeing with iHaubois et al. 1 (12009). The combination of 



this angular diameter and bolometric flux and distance yields an effective temperature 
Teff = 3590 ± 55 K and 7? = 955 ± 217 R^, respectively. 

The next step is to fit the limb-darkening law Eq. 1 to the interferometry data. 
We try to determine limb-darkening coefficients from the first two lobes of the visibil- 
ity curve, but the resulting limb-darkening fit predicts large intensities near the limb, 
7(0) /2'K > 0.7, which is greater than one-half of the central intensity. This suggests 
that the second lobe of the visibility curve also contains light from circumstellar excess. 
Therefore, we attempt to fit the limb-darkening law to the first lobe only, even though 
the first lobe of the visibility curve is a weak measure of the intensity profile. We also 
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Figure 3. Con-elation between limb-darkening coefficients A and B from 
spherically-symmetric model stellar atmospheres with effective temperatures Teif = 
3400 - 3800 K. 



consider the predicted limb-darkening coefficients for models with effective tempera- 
tures ranging from 3400 K to 3800 K. In Fig.|3j we plot the limb-darkening coefficients 
from model atmospheres in that temperature range. In the range of A = - 3 there is a 
tight correlation between the coefficients of the form B = -1.27A -i- 1.50 . Using this 
correlation, we reduce the number of degrees-of-freedom. 

When we fit the observed first lobe of the visibility curve for the limb-darkening 
coefficients, we find A = 1.36 ± 0.33, which leads to B = -0.22 ± 0.43. Using these 
limb-darkening coefficients with the correlations, Eq. 2, and 3, where = 4.338 x 
lO"'^, Dp - 0.73704, C/ = 0.01966, and D/ = 1.0601, we find a value of RJM^ = 
82.17+}[^2;jo/Mo for Betelgeuse, hence for/? = 955±217 /?o we get M - 11.6!^° Mq. 
Using the mass and radius, we predict the gravity to be logg = -0.45^° }^. 

As a check, we compare these results to limb-darkening fits to model atmospheres, 
the "nearest" model to the predicted effective temperature, gravity and mass is a T^g = 
3600 K, logg = -0.5, and M = 10 M© model atmosphere. Its limb-darkening coeffi- 
cients are A = 1.63 and B = -0.55 with a ratio RJM^ ^ 93.32 Rq/Mq. This ratio of 
Rt /Mh. agrees with the predicted ratio for Betelgeuse. 
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6. Summary 

In this work, we use a flux-conserving Umb-darkening law, sphericany-symmetric model 
steflar atmospheres computed with the SAtlas code, and interferometric observations 
to probe the fundamental parameters of Betelgeuse. This novel method computed one 
of the first measurements of the mass of Betelg euse. The predicte d mass is less than 



that determined by stellar evo lution calculations (jPolan et al.ll2008r) . but it is consistent 



with the gravity suggested by iLobel & Dupreel ( 2000 ). 



The predicted mass is still very uncertain, the two main sources of uncertainties 
are the distance to Betelgeuse and the circumstellar flux excess apparent in the interfer- 
ometric observations. In spi te of the large uncertainties, the mass prediction is precise 



enough to agree with lSmith et al. (2009) that Betelgeuse will explode as a Type IIP su- 
pernova, and may even act as a constraint for stellar evolution models and input physics, 
such as convective core overshooting. 

It is also worth noting that interferometry at shorter wavelengths may be a more 
effective constraint for our method. At shorter wavelengths, interferometry has greater 
angular resolution, as well as avoiding circumstellar emission in the infrared. 
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